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Abstract 

In the genome of the biotrophic plant pathogen Ustilago maydis, many of the genes coding for secreted protein effectors 
modulating virulence are arranged in gene clusters. The vast majority of these genes encode novel proteins whose 
expression is coupled to plant colonization. The largest of these gene clusters, cluster 19A, encodes 24 secreted effectors. 
Deletion of the entire cluster results in severe attenuation of virulence. Here we present the functional analysis of this 
genomic region. We show that a 19A deletion mutant behaves like an endophyte, i.e. is still able to colonize plants and 
complete the infection cycle. However, tumors, the most conspicuous symptoms of maize smut disease, are only rarely 
formed and fungal biomass in infected tissue is significantly reduced. The generation and analysis of strains carrying sub- 
deletions identified several genes significantly contributing to tumor formation after seedling infection. Another of the 
effectors could be linked specifically to anthocyanin induction in the infected tissue. As the individual contributions of these 
genes to tumor formation were small, we studied the response of maize plants to the whole cluster mutant as well as to 
several individual mutants by array analysis. This revealed distinct plant responses, demonstrating that the respective 
effectors have discrete plant targets. We propose that the analysis of plant responses to effector mutant strains that lack a 
strong virulence phenotype may be a general way to visualize differences in effector function. 
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Introduction 

U. maydis is a biotrophic fungal pathogen causing smut disease in 
maize. To cause disease, haploid cells of compatible mating type 
need to fuse on the plant surface and develop an infectious 
dikaryon [1,2]. Upon perception of appropriate surface cues [3], 
the dikaryon differentiates non-melanized appressoria that pene- 
trate plant cells direcdy, presumably aided by local secretion of 
lytic enzymes [4] . During penetration, the host plasma membrane 
invaginates and encases the fungal hyphae, a feature typical for 
biotrophs. This establishes an extended interaction interface for 
the exchange of signals and nutrients [5,6]. Initial intracellular 
biotrophic growth of U. maydis is followed by intercellular growth 
during later stages of the infection, concomitant with massive 
proliferation in mesodermal tissue close to the veins. At this 
developmental stage, huge fungal aggregates form in cavities 
between plant cells followed by differentiation of ornamented 
diploid spores [6]. Fungal proliferation coincides with plant cell 
enlargement and resumption of mitotic divisions [6] . U. maydis can 
infect and cause symptoms on all above ground maize organs, with 
the infection staying locally confined. This is in contrast to related 
smut fungi that show systemic spread throughout the plant but 
produce symptoms only in the male and female inflorescences 
[7,8]. 



During initial contact of U. maydis with the plant leaf and 
presumably triggered by fungal pathogen-associated molecular 
patterns (PAMPs), a number of plant defense genes are induced. 
This upregulation disappears during fungal penetration, suggest- 
ing that these initial defense responses are actively suppressed by 
the fungus during the plant colonization stages [9]. Also, several 
genes associated with suppression of plant cell death are induced. 
One of these, the maize cystatin CC9, has recendy been 
functionally analyzed [10]. Silencing of CC9 enhanced maize 
defense gene expression and upon infection with U. maydis a 
hypersensitive response was observed. CC9 was shown to suppress 
apoplastic cysteine protease activity, illustrating that CC9 is a 
novel compatibility factor for the biotrophic interaction of maize 
with U. maydis [10]. After plant colonization the most dramatic 
transcriptional changes in the host affect hormone signaling, 
induction of antioxidants, secondary metabolism, as well as a block 
in the transition from a juvenile sink tissue to a mature, 
photosynthetically active source tissue normally observed during 
leaf establishment [9]. The latter supports the observation that U. 
maydis is able to colonize young meristematic maize tissue, but is 
unable to infect differentiated source tissue [11]. The changes in 
plant gene expression observed after host colonization are likely to 
be brought about by secreted fungal effector molecules. The 
genome of U. maydis encodes about 300 novel secreted effectors 
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Author Summary 

In this study, we provide the first step to the functional 
analysis of the largest gene cluster in the Ustilago maydis 
genome encoding 24 secreted effectors. While the 
deletion of the entire cluster dramatically affected tumor 
formation and abolished anthocyanin induction, only one 
of the genes had a large contribution to tumor formation, 
while another effector gene was primarily responsible for 
the anthocyanin induction. Unexpectedly, the cluster 
mutant could still colonize plants and complete the life 
cycle, i.e. behaves like an endophyte. Despite only small 
contributions to tumor formation, individual effector 
mutants caused distinct plant responses, suggesting that 
they affect discrete plant processes. On these grounds we 
are proposing to use plant responses as a general readout 
to assess and compare effector gene function. 

that are upregulated during plant colonization and largely lack 
known InterPro domains. Of these a significant percentage is 
arranged in gene clusters and deletion of entire clusters can have 
dramatic effects on virulence [12,13]. Genome comparisons with 
the related smut fungi Sporisorium reilianum and U. hordei revealed 
that the majority of secreted effectors also exist in these relatives. 
With respect to conservation, effectors fall in two classes: 
approximately 34% are highly conserved in all smut fungi 
sequenced so far [8] and the remainder are poorly conserved, 
reflecting the arms race with the host. Furthermore, in U. hordei the 
tight clustering of effector genes seen in U. maydis and S. reilianum is 
largely disrupted [7,8]. Effector genes in U. maydis are plant- 
induced [12], and work by Skibbe et al. [14] has revealed that the 
expression of some effector genes is tissue specific, i.e. is different 
when U. maydis colonizes seedlings, adult leaves or tassel, and the 
need for the corresponding effectors may be restricted to the 
respective tissues. So far the function of only very few of the many 
novel effectors has been elucidated in the U. maydis/ 'maize 
pathosystem. Pepl, a conserved effector is needed for penetration 
[15] and affects plant defense responses by inhibiting apoplastic 
plant peroxidases [16]. Pit2, another conserved effector affects host 
defense responses [17] through inhibition of apoplastic plant 
cysteine proteases [18]. Cmul is a secreted chorismate mutase that 
is taken up by plant cells and lowers salicylic acid (SA) levels in 
infected tissue through metabolic priming [19]. 

Here we describe the beginning of the functional analysis of 
cluster 19A, the largest effector gene cluster in U. maydis. In a 
previous study, the entire cluster 19A comprising 23 putative 
effector genes in a 40 kb genomic region was deleted [12]. Cluster 
1 9A mutants were severely attenuated in virulence and except for 
some ligula swelling, the mutants rarely elicited tumor formation 
[12]. In this study, we map the most important effectors for 
seedling infection in this cluster, and show that most strains deleted 
for individual effector genes show only minor reductions in 
virulence but elicit distinct plant responses. 

Results 

Mapping of cluster 19A genes contributing to tumor 
formation and anthocyanin induction after seedling 
infection 

Cluster 19A was originally predicted to encode 23 secreted 
effectors [12]. The manual reannotation based on comparison 
with S. reilianum and U. hordei now predicts the presence of 24 
effector genes, plus one gene related to a reverse transcriptase 
(um05313) presumably originating from a retrotransposon and one 



pseudogene (um05316) (http://mips.helmholtz-muenchen.de/ 
genre/proj/ustilago) (Figure 1A). 

Among the effector genes in cluster 19A, we detect five gene 
families based on amino acid sequence similarity (Figure 1A, 
Figure 2 and Figure SI). Um05294, Um05295, Uml2302, 
Uml0553 and Uml0554 display between 34-48% similarity at 
the amino acid sequence level (Figure 1A). Another family 
(Figure 1A) comprises genes um05299, um05300 and um05301. 
The respective proteins show between 34 — 4 1 % similarity. The two 
effector proteins encoded by the adjacent genes um05305 and 
um(>5306 display 30% amino acid similarity (Figure 1A). A three 
gene family codes for Um05309, Um05310, Um05311 with 35- 
86% amino acid sequence similarity and the largest family is 
comprised of Um05312, Um05314, Uml0557, Um05317, 
Um05318 and Um05319 with 33-51% amino acid similarity 
(Figure 1A). Outside of cluster 19A the U. maydis genome does not 
contain paralogs to any of these gene families. The 24 effectors 
encoded by cluster 19A do not contain recognizable protein 
domains nor do they display a characteristic spacing of cysteine 
residues described for several other U. maydis effectors [13]. 
However, orthologs for most of these genes are found in the 
genomes of S. reilianum and U. hordei [7,8] (http://mips.helmholtz- 
muenchen.de/ genre/ proj/ sporisorium/) (http:/ / mips.helmholtz- 
muenchen.de/genre/proj/MUHDB/) (Figure 2). Published ex- 
pression data for cluster 19A genes from different U. maydis 
infected tissues [12,14] are compiled in Figure 2. These studies 
revealed that except for two genes where expression could not be 
detected, genes in cluster 19A are differentially induced when 
different plant organs are colonized (Figure 2). Furthermore, only 
three of the cluster 1 9 A genes are downregulated when the central 
regulator for pathogenic development, the bE/bW complex, is 
switched off during biotrophic development [20] (Figure 2). This 
illustrates, that the individual cluster genes are plant induced but 
do not appear to be co-regulated. 

Given the strong virulence phenotype of cluster 19A mutants 
[12], we examined whether the cluster 19A mutant is impaired in 
biotrophic growth. To this end, we compared biotrophic 
development of the mutant strain SG200A19A and the progenitor 
strain SG200. To our surprise, SG200A19A formed appressoria 
on the plant surface, proliferated inside plant tissue and 
differentiated teliospores, at a late time point comparable to 
infections with SG200. This explains, why spores had not been 
detected in the previous study [12] and illustrated that the 19A 
mutant could undergo biotrophic development and complete the 
life cycle (Figure IB). With respect to macroscopic symptoms, 
tumor formation was abolished and instead chlorotic and necrotic 
areas became apparent (Figure 1C). In addition, anthocyanin 
accumulation was completely abolished in mutant infected tissue, 
while prominent anthocyanin stained regions were observed in 
leaves infected by SG200 (Figure 1C). 

To identify the genes contributing to the phenotype of the 
cluster 19A mutant, a series of strains was generated carrying sub- 
deletions of cluster 19A. At first, we divided cluster 19A into two 
parts which we designate 19A-1 (left half region of cluster 19A) and 
19A-2 (right half region of cluster 19A) (Figure 1A) and generated 
respective deletion mutants. When tested for phenotype after 
seedling infection in comparison to SG200, SG200A19A-1 showed 
a dramatic reduction of tumor formation and loss of anthocyanin 
accumulation. The effects were comparable to infections with 
SG200A19A. Conversely, SG200A19A-2 was only weekly atten- 
uated in virulence and was able to elicit anthocyanin accumulation 
(Figure 1C and Figure 3). To determine possible differences in the 
efficiency of plant colonization by SG200 and the different derived 
mutant strains, we analyzed fungal biomass in colonized tissue by 
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Figure 1. Phenotypic analysis of cluster 19A mutants. A) Schematic representation of the cluster 19A region in the genome of U. maydis. 
Arrows indicate open reading frames (ORFs), labeled by the respective gene numbers. Related gene families are indicated in orange, yellow, blue, 
green and brown and the percentages of amino acid similarity are given below. Unique effector genes are shown in grey. Genes encoding proteins 
with a predicted secretion signal are depicted in orange, yellow, blue, green, brown or grey while open arrows indicate genes encoding proteins 
without a prediction for secretion. The left half of cluster 19A was designated as 19A-1 and the right half was designated 19A-2 as indicated. B) 
Biotrophic development of SG200 and SG200A19A. Maize seedlings were infected with SG200 or SG200A19A and development was monitored 
microscopically. At 1 dpi penetrating hyphae on the leaf surface were visualized by calcofluor staining. Proliferating biotrophic hyphae were 
observed after chlorazole-black E staining at 13 dpi. At 30 dpi mature teliospores were observed indicated by yellow arrows. Bar=100 \im. C) 
Macroscopic symptoms of representative leaves infected with SG200, SG200A19A, SG200A19A-1 and SG200A19A-2 at 12 dpi. Note the absence of 
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anthocyanin pigment accumulation in infections with SG200A19A and SG200A19A-1. D) Quantification of fungal biomass by qPCR. Genomic DNA 
was extracted from the leaves infected with SG200, SG200A19A, SG200A19A-1 and SG200A19A-2 at 2 and 6 dpi and used for qPCR. Relative fungal 
biomass was calculated by the comparison between U. maydis peptidylprolyl isomerase gene (ppi) and Z. mays glyceraldehyde 3-phosphate 
dehydrogenase gene {GAPDH). Error bars indicate standard deviation. ** p<0.01. p-values were calculated by Student's f-test. 
doi:1 0.1 371 /journal.ppat.1 003866.g001 



quantitative real time PCR using U. maydis peptidylprolyl 
isomerase (ppi) and Z- ma y s glyceraldehyde 3-phosphate dehydro- 
genase [GAPDH] as reference genes for quantifying fungal and 
plant biomass, respectively [15,21]. At 2 days post infection (dpi), 
fungal biomass of SG200A19A, SG200A19A-1 and SG200A19A- 
2 was comparable to the SG200 infection (Figure ID, left panel). 
At 6 dpi, however, fungal biomass of SG200A19A and 
SG200A19A-1 was lower compared to SG200 (Figure ID, right 
panel), suggesting growth defects of these mutant strains at this 
later time point where massive fungal proliferation inside tumors is 
observed in SG200 infected tissue. Fungal biomass of 
SG200A19A-2 lacking the right half of cluster 19A was not 
significandy different from SG200 (Figure ID), consistent with the 
weak effect of this deletion on virulence. These results indicated 
that the major effector genes responsible for the phenotype of the 
cluster 19A mutant after seedling infection must reside in the left 
half of cluster 19A. 

To determine the contribution to virulence of the 14 genes 
located in the left half of cluster 19A, we divided the region into 
four parts, 19A-la, 19A-lb, 19A-lc and 19A-ld making sure that 
existing gene families were deleted simultaneously (Figure 3A). 
Respective deletion mutants were generated and assayed for tumor 
formation and ability to elicit anthocyanin accumulation 
(Figure 3B,C). Of these four mutants SG200A19A-la and 
SG200Al9A-lc did show a small reduction of tumor formation, 
but this was not statistically significant (Figure 3B). The deletion of 
19A-lb and 19A-ld significandy lowered tumor formation with 
19A-ld showing the strongest effect (Figure 3B). Anthocyanin 
accumulation was abolished by the deletion of 19A-lc and was 
unaffected in infections with the other sub-deletion mutants 
(Figure 3C). A triple deletion generated by combining 19A-lb, 
19A-lc and 19A-ld deletions resulted in a mutant strain 
(SG200A19A-lbcd) that was severely reduced in tumor formation 
and failed to induce anthocyanin, comparable to the deletion 
strain lacking the left half of cluster 19A (SG200A19A-1). This 
firmly establishes that region 19A-la does not measurably 
contribute to virulence (Figure 3B). 

Based on the finding that the um05318 gene in U. maydis cluster 
19A has three paralogous genes at a syntenic position in cluster 
19A of S. reilianum [7] (Figure 2) that contribute to virulence in S. 
reilianum (H. Ghareeb and J. Schirawski, personal communication), 
we decided to delete the two rightmost genes in U. maydis cluster 
19A (19A-2e, Figure 3A). This was done also to rule out that the 
lack of a virulence phenotype when deleting the entire right half 
(A19A-2, Figure 3B) is caused by balancing positive and negative 
effects of effector gene deletions on virulence [12]. The resulting 
strain designated SG200Al9A-2e showed weakly attenuated 
tumor formation that was comparable to strain SG200A19A-2 
deleted for the entire right half of the cluster and this was again not 
statistically relevant (Figure 3B). 

Identification of major effector genes contributing to 
cluster 19A deletion phenotypes 

To identify the major effector genes that contribute most 
strongly to the virulence phenotype of sub-clusters 19A-lb, 19A- 
lc, 19A-ld and 19A-2e, we initially generated overlapping sub- 
deletions and tested them for virulence (not shown). For example, 
the cluster 19A-ld region was subdivided into a double deletion of 



um05305/um05306 and a double deletion of um05306/uml0556, 
respectively. These double mutants were tested for virulence and 
in this case only the double deletion of um()5306 Vuml 0556 was 
affected in virulence, allowing the conclusion that uml 0556 is the 
responsible gene (not shown). This was then followed up by single 
deletions of the genes identified in this approach (Figure 4A). With 
respect to 19A-lb comprising five related genes designated tinl-1 to 
tinl-5, we were unable to identify the respective major individual 
effector gene(s), as the observed effects of further sub-deletions on 
tumor formation could no longer be assessed as being statistically 
relevant (data not shown). This suggested that these five members 
of the same gene family (Figure 2) contribute weakly but additively 
to virulence. Complementation experiments where all 5 genes 
were re-introduced in strain SG200A19A-lb revealed that the 
weak virulence phenotype of the 19A-lb deletion could be 
complemented (Figure 4B). For the region 19A-lc the deletion of 
um05302 (designated tin2) showed a comparable reduction in 
virulence to the deletion of the entire 19A-lc region. In addition, 
the deletion oitin2 abolished anthocyanin induction (Figure 4B,C). 
The introduction of a single copy of tin2 into SG200Atin2 
complemented the weakly reduced tumor formation as well as 
anthocyanin accumulation (Figure 4B,C). The single deletion of 
uml0556 (designated tinS) significantly affected tumor formation 
(Figure 4B). Also in this case, tumor formation could be restored to 
a level comparable to SG200 by introducing tin3 in single copy 
(Figure 4B). After infection, the single gene deletion mutant of 
um05318 (designated tin4) showed a lower incidence of plants with 
stronger disease symptoms than SG200 infected plants (Figure 4B). 
Upon complementation, this disease category was increased 
compared to SG200Atin4 and more plants showed heavy tumors 
and were stunted or dead (Figure 4B). This suggests a weak 
contribution of tin4 to tumor formation. With respect to um05319 
(designated tin5), the single gene deletion had minor effects on 
virulence and the reintroduction of the gene did not significandy 
change the disease scores (Figure 4B). 

Plant responses to cluster 19A mutants 

The analysis of single effector gene mutants for cluster 19A 
revealed in general, that deletions of individual genes (with the 
exception of tinS) had only minor or statistically non-substantial 
effects on virulence (Figure 4B), suggesting that the strong 
virulence defect observed in the entire cluster 19A deletion is 
due to additive effects and/ or concerted action of several effectors. 
To visualize this and to obtain evidence whether individual 
effectors target distinct plant processes, we decided to analyze the 
plant responses to infection by the 1 9A deletion strain as well as to 
several single effector gene mutants on the transcriptome level. 
Maize seedlings were infected by SG200A19A, SG200A19A-lb, 
SG200Atin3, SG200Atin4 and SG200Atin5. RNA was extracted 
from infected plant material harvested at 4 dpi, a time point where 
the individual mutants should not differ in fungal biomass as 
assessed from the analysis of sub-deletions (Figure ID). Three 
biological replicates were prepared and analyzed by Affymetrix 
maize genome microarrays. For technical reasons, the expression 
data for plants infected with the SG200Atin2 mutant strain could 
not be included in this comparative transcriptome analysis. 

Maize gene expression profiles of tissue infected with cluster 
mutant strains were compared to profiles of SG200 infected and 
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Figure 2. Expression data compilation for cluster 19A genes and presence of orthologous genes in other smut fungi. Expression data 
for cluster 19A genes were compiled from the following sources: a) Gene expressed in young seedlings [14], b) Gene expressed in tassel [14], c) Gene 
expressed in tumor tissue [12], d) Gene down-regulated in response to b inactivation during biotrophic development [20] and nd) no expression 
detected. Orthologs in S. reilianum and U. hordei are listed, and in these columns ND indicates that no orthologs were detected at a cut-off of e-value 
1 0e~ 10 . The adopted color scheme corresponds to the scheme used in Figures 1 , 3 and 4. Please note that the order of genes was changed in several 
cases and does not follow the order in cluster 19A (*). 
doi:1 0.1 371 /journal.ppat.1 003866.g002 
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mock-treated plants, which had been generated in our previous 
study on the transcriptional responses of maize to U. maydis and 
had been used as reference in the analysis of plant responses to 
pepl and pit2 effector mutants [9,15,17]. RMA-normalized 
microarray data were then subjected to a one-way ANOVA and 
contrast gene lists were generated using a fold change of ±2 and a 
corrected /j-value of 0.05 as cutoffs (Table SI). Expression of a set 
of 13 genes differentially regulated after infection with different 
mutant strains was subsequently analyzed by quantitative real-time 
PGR (qRT-PCR) with RNA from independently generated 
infected plant material, and this allowed validating the array 
results (Figure S2). 

Compared to SG200 infected samples, 1816 maize genes were 
differentially regulated in response to SG200A19A (Table SI). A 
hierarchical clustering of these 1816 genes was performed for the 
whole data set to visualize the relations between the transcriptional 
responses to the individual U. maydis strains (Figure 5A). As 
expected, the maximal distance in gene expression was between 
SG200 infections and infections by the SG200A19A mutant, 
which caused only very weak disease symptoms and thus displayed 
the highest similarity to the mock-inoculated plants (Figure 5A). 
On the other hand SG200Al9A-lb infections showed highest 
similarity to SG200 infections, illustrating that the 5 tinl effector 
genes have only a weak contribution to plant responses, which is in 
line with their weak effect on virulence. Profiles of plant responses 
to strains carrying tin gene deletions clearly discriminated 
SG200Atin3 and SG200A19A-lb, while responses to SG200Atfn4 
and SG200Atin5 infections were not separated by the hierarchical 
clustering (Figure 5A), indicating that similar responses were 
elicited by tin4 and tin5 mutants. 

We also combined the transcriptional responses elicited by each 
of the four individual tin mutant strains to be able to compare this 
to the response elicited by SG200A19A, the strain carrying the full 
deletion of cluster 19A and to reveal contributions of genes not 
deleted individually. In total, 1513 maize genes were differentially 
regulated by the four tin mutants compared to SG200 infections, 
while 1816 genes were differentially regulated after infection with 
the cluster deletion mutant SG200A19A (Figure 5B). Interestingly, 
a comparison of these combined "responses to individual tin gene 
deletions" to the SG200A19A responsive genes showed only a 
partial overlap of differentially regulated transcripts (Figure 5B). 
Of the 726 genes differentially expressed in response to the whole 
cluster mutant but not detected in the "responses to individual tin 
gene deletions" 352 genes were induced and 374 genes were 
repressed (Figure 5B). Among the induced transcripts, particularly 
biotin synthesis genes were induced specifically in SG200A19A 
infected tissue, while plant cellulose synthesis genes were 
downregulated after infection by the A19A deletion strain (Table 
S2 and Figures S3 and S4). In addition, several anthocyanin 
biosynthesis related genes were downregulated after infection with 
the cluster 19A mutant while they were not included as 
differentially expressed in "responses to individual tin gene 
deletions" (Table S2). This most likely indicates a contribution 
of the Tin2 effector to anthocyanin induction, and reflects that 
plant responses to the tin2 deletion strain could not be included in 
"responses to individual tin gene deletions" for technical reasons. 
Of the 1090 "shared differentially regulated plant genes", 288 
genes were induced, while 802 genes were downregulated 
(Figure 5B). Induced genes comprised pathogen response genes 
such as PR4, PR5 and several oxidases, demonstrating an elevated 
plant defense in response to tin gene and cluster 19A mutant 
strains (Table S2). Downregulated transcripts were strongly 
enriched for genes involved in DNA-metabolism and DNA- 
modification, particularly histones and DNA-methyltransferases 



(Table S2 and Figure S5). This most likely reflects the reduced tumor 
formation observed in all the mutants compared to SG200. On the 
other hand, 423 genes were differentially regulated after infections 
with the tin mutant strains (Figure 5B) but these were not differentially 
regulated after infection with the cluster 19A mutant. Among these 
423 genes, several chitinases and peroxidases were found (Table S2). 

To get clues on the possible roles of the tin genes during host 
colonization we next visualized the differentially regulated plant 
genes in response to individual tin mutants. This analysis revealed 
that 476 maize genes were commonly regulated by all four 
mutants, while 1027 genes only responded to a subset of tin 
mutants (Figure 5C). The smallest number of specifically regulated 
maize genes was found after infections with the SG200A19A-lb 
mutant, which lacks the five related tinl effector genes, i.e. this 
mutant shared 73% (476 out of 65 1) of differentially regulated genes 
with all other strains. Amongst the 66 plant genes that were 
specifically regulated after SG200A19A-lb infections, four maize 
endochitinase genes were significantly induced compared to SG200 
infected leaves. In addition, transcript levels of two salicylic acid 
binding proteins and peroxidase- 12, which was found to be involved 
in the maize apoplastic oxidative burst [16], were induced (Table 
S3). This suggests that tinl genes modulate basal defenses. 

The tin3 mutant specifically affected the differential regulation of 
70 maize genes (Table S3). Conspicuously, sucrose synthase and 
several transcription factors including auxin-response factors were 
induced suggesting a link to the reduced ability for tumor formation. 

After infections by tin4 and tin5 mutant strains the majority of 
differentially expressed maize genes were shared (Figure 5C), which 
is in accordance with the hierarchical clustering result which places 
these two strains closely together (Figure 5A). Nevertheless, 167 
genes were only differentially regulated by the tin4 mutant and 1 10 
maize transcripts were differentially regulated in infections by the 
tin5 mutant (Figure 5C and Table S3). Although tin4 and tin5 
mutants displayed an indistinguishable virulence behavior, they 
elicited distinct molecular responses in maize. This appears to be a 
general feature and for example also holds for a strain carrying the 
19A-lb deletion (Figure 3B) where virulence is only weakly 
attenuated. When the 10 most strongly induced host genes are 
compared after SG200A19A-lb infection, only one gene was 
among the top 10 genes upregulated after infection with 
SG200Atin3, SG200Atin4 or SG200Atin5, respectively (Table SI). 

Gene ontology enrichment analyses were performed for the plant 
genes upregulated after infection by the three single tin gene deletion 
strains SG200Atin3, SG200Atin4 and SG200Atin5 (Figures S6, S7, 
S8). Differentially enriched functions for the individual effector 
deletion strains were visualized and in addition genes corresponding 
to the enriched functions are listed (Figure 6 and Table S4). This 
revealed that statistically distinct processes were induced in plants 
infected by the different mutants. For SG200Atin3, oxidoreductases 
and carbonate dehydratases were enriched, while SG200Atin4 
infection induced plant genes significantly enriched for functions 
involved in iron ion binding. In SG200Atin5 infection, on the other 
hand, lipoxygenases and serine-carboxypeptidases were induced, 
which were not as highly upregulated in the other mutants (Figure 6 
and Table S4). These results illustrate that plant responses to 
effector mutants which show no or only weak reductions in 
macroscopic symptoms can be highly specific and can be used to 
describe and discriminate mutant phenotypes. 

Discussion 

In this communication, we have dissected the largest U. maydis 
effector gene cluster 19A, identified the most relevant effectors for 
seedling infection and present evidence, that individual effectors 
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Figure 3. Mapping of effector genes responsible for the virulence phenotype of cluster 19A mutant. A) Schematic representation of 
analyzed sub-deletions. The region of 19A-1 was divided into 4 parts, 19A-1b, 19A-1a, 19A-1c and 19A-1d according to the distribution of gene 
families. In the second half of the cluster (1 9A-2) only the 1 9A-2e region was deleted. B) Quantitative evaluation of virulence of cluster 1 9A mutants. 
The scoring scheme follows the one described by Kaemper et al. [12] with severity of virulence corresponding to the color intensity. The numbers on 
top of respective bars indicate the total number of infected plants. Results from at least 3 independent experiments are combined. Dots above 
respective panels indicate that at least one of symptoms (ligula swelling, small tumors, normal tumors, heavy tumors, stunted or dead) was 
significantly changed in the mutant relative to SG200 by Student's f-test. C) Quantification of infected plants showing anthocyanin accumulation. In 
the same experiments as in B) the percentage of plants showing anthocyanin accumulation was scored. 
doi:1 0.1 371 /journal.ppat.1 003866.g003 



target distinct processes in the host plant. It was published before, 
that the deletion of cluster 19A abolishes tumor formation [12]. 
We now show that this dramatic phenotype is not associated with a 
block in biotrophic development. The cluster 1 9A mutant was still 
able to complete the life cycle up to the formation of teliospores. 
However, massive fungal proliferation observed at later stages in 
tumor tissue in infections with wild type strains was absent, 
suggesting that the effectors in this cluster are responsible for 
tumor induction either direcdy or indirectly but not for growth per 
se in the infected tissue. The analysis of plant responses elicited by 
the whole cluster mutant (SG200A19A) in comparison to 
responses to the progenitor strain SG200 revealed that 1816 of 
the 13,339 maize genes represented on the chip were differentially 
regulated. The analysis of the maize transcriptome changes 
observed for the individual effector mutant infections revealed 
that about 60% of these changes were shared by all mutants, 
suggesting that they are unspecific. Another aspect of the 
comparative transcriptome profiling was the finding that there is 
incomplete overlap between the genes affected by the individual 
mutants and genes altered in their regulation when the entire 
cluster is deleted. This could indicate that the effects of individual 
effector deletions cease to be visible when the entire cluster is 
deleted, i.e. the dramatic phenotype of the cluster deletion might 
bury the more subtle physiological changes caused by individual 
mutants. Our study has also revealed that effectors not studied 
individually because of their undetectable contribution to 
virulence after seedling infection, may profoundly affect the 
metabolic activity of the infected tissue. An example is the 



observed upregulation of biotin biosynthesis after infection with 
SG200A19A which is not observed in any of the tin mutant 
infections. In Arabidopsis thaliana it has been shown that biotin is 
critical for suppressing spontaneous cell death [22]. The fact that 
biotin biosynthesis appears upregulated may indicate a direct 
involvement of a specific effector for maintaining a certain level of 
this essential cofactor. Alternatively, biotin upregulation could be a 
secondary effect that allows the cluster 1 9A mutant to grow in the 
infected tissue. Future array analyses with mutants of cluster 19A 
not studied here with respect to the plant responses they elicit 
should allow to separate primary and secondary effects and allow 
to uncover the effector responsible for the regulation of glycolysis 
and biotin biosynthesis, respectively. 

Among the common genes differentially expressed after 
infection with all tin mutant strains, we observe an enrichment 
of upregulated plant defense genes and downregulation of genes 
involved in DNA metabolism. This is likely to reflect insufficient 
suppression of plant defenses due to reduced fungal proliferation 
(or the absence of certain effectors) and reduced plant tumor 
formation, respectively. Similarly, genes for photosynthesis com- 
ponents were consistently higher expressed after infections with all 
mutant strains compared to SG200 infections. This is unlikely to 
reflect an induction of photosynthesis during mutant infections but 
presumably results from an incomplete shutdown of photosynthe- 
sis, usually observed during U. maydis wild type infection [9] . The 
reduction in plant cell wall biosynthesis gene expression after 
infection with the cluster 19A mutant, which is not seen in 
infections with individual tin mutants, likely reflects the reduced 
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Figure 4. Virulence of fmgene deletion mutants and complemented strains. A) Schematic representation of the deleted tin genes in cluster 
19A. tin1-1 to tinl-5 correspond to the 5 genes represented by 19A-1 b, the other tin genes lie in region 19A-1c (tin2), region 19A-1d (tin3), and region 
19A-2e (tin4 and tinS). B) Quantitative evaluation of virulence of tin gene deletion mutants and complemented strains. Virulence of the displayed 
deletion mutants and respective complementation strains (indicated as A and C, respectively below each bar) was determined as described in 
Figure 3B. C) Quantification of infected plants showing anthocyanin accumulation after infection with tin2 mutant and control strains. Plants were 
infected with the indicated strains and the percentage of plants showing anthocyanin accumulation was scored 12 dpi. 
doi:10.1371/journal.ppat.1003866.g004 



ability of the cluster mutant to induce tumors containing enlarged 
plant cells [15], while individual tin mutants can still induce 
tumors. 

With a minor impact on tumor formation, Tin2 was specifically 
responsible for anthocyanin accumulation in infected tissue while 
all other individual mutants showed anthocyanin induction. The 
transcriptome analysis revealed that several anthocyanin biosyn- 
thesis genes were upregulated in all the individual tin mutant 
strains but not after infection with SG200A19A lacking the whole 
cluster including tin2 (Table SI). This suggests that Tin2 is directly 
responsible for inducing these anthocyanin biosynthesis genes. 
Anthocyanin has been hypothesized to have a protective role 
against abiotic stresses [23], and can be induced after biotic stress, 
although its role here is unclear [24]. 

While the virulence assays for tin4 and tin5 mutant strains were 
largely uninformative because of limited assay sensitivity, the 
transcriptome analysis revealed that tin4 and tin5 mutants elicited a 
series of plant responses that were mutant specific, but in addition 
318 differentially regulated plant genes were differentially regu- 
lated by tin4 as well as tin5 mutant strains. Based on the fact that 
Tin4 and Tin5 share 19% identity and 39% amino acid similarity 
this could indicate that these effectors are in the process of 
diversification to different functions (167 tin4 specifically regulated 
transcripts, 1 1 0 tin5 specifically regulated transcripts) while still 
maintaining some of the original common functions (318 
commonly regulated transcripts). This interpretation would also 



make sense in view of the fact that um05312, um05314, uml0557 
and um05317 which are also related to tin4 and tin.5 (Fig. S1E) were 
not individually deleted in this study. If all of the genes in this 
family had redundant functions, we would not have expected to 
see differences in the host responses to the tin4 and tin5 mutants. 
The gene ontology enrichment analysis showed that several genes 
involved in iron metabolism/uptake were specifically upregulated 
when tin4 is missing. Elevated iron availability may directiy affect 
the activity of the respiratory burst oxidase requiring a heme 
prosthetic group to generate superoxide [25]. Such effects on plant 
defense have also been described after infections by Erwinia 
crysanthani [26] and Blumeria graminis f. sp. tritici [27]. After tin5 
mutant infections lipoxygenases and serine-carboxypeptidases 
were specifically upregulated, and both types of enzymes have 
been associated with defense [28,29]. OsBISCPLl, a serine- 
carboxypeptidase from rice, was up-regulated in incompatible 
interactions between rice and the blast fungus, and was implicated 
in regulation of defense responses from heterologous expression 
studies [30]. In S. reilianum three orthologs of tin4 are present 
[sr 10075, sr 10077 and sr 10079; Figure 2). The simultaneous 
deletion of the neighboring genes srl0073, srl0075, srl0077 and 
sr 10079 weakly affected S. reilianum virulence (H. Ghareeb and J. 
Schirawski, personal communication), similar to what we observe 
for the tin4, tin5 double mutant of U. maydis. As in U. maydis, the left 
half of cluster 19A contributes most strongly to virulence in S. 
reilianum (H. Ghareeb and J. Schirawski, personal communication). 
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Figure 5. Plant responses to cluster 1 9A mutant strains. A) Hierarchical clustering visualizing the relative expression of the 1816 maize genes 
transcriptionally regulated 4 days after infection by U. maydis strain SG200A19A. X-axis depicts clustering of the microarray samples, with three 
biological replicates per strain. Y-axis shows clustering of the regulated maize transcripts based on the similarity of their expression patterns. B) Venn 
diagram showing numbers of maize transcripts differentially expressed 4 days after infection with SG200A19A (red circle) as compared to the 
combined deregulated transcripts of all individual tin mutants (green circle). For the tin mutants, all genes regulated in at least one of the mutant 
infections have been combined. Red numbers indicate induced expression; numbers of down-regulated transcripts are given in blue. C) Venn 
diagram showing maize transcripts differentially expressed 4 days after infection by the U. maydis mutants SG200A19A-1b (A19A-1b), SG200Atin3 
(Atin3), SG200Atin4 (Atin4), SG200Atin5 (Atin5), respectively. 
doi:1 0.1 371 /journal.ppat.1 003866.g005 



Based on the observation that cluster 19A effector genes of U. 
maydis are not essential for tumor formation in tassel (although 
tumor size in tassel was reduced after infection with SG200A19A) 
[14], we consider organ-specificity of effector function more likely 



to explain this finding. As S. reilianum does not induce tumors in 
leaves and develops disease symptoms only in the cob and in the 
tassel, this species may not need effector genes like tin3 for tumor 
induction in vegetative tissues of the maize plant. Consistent with 
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this is our observation that the U. maydis cluster 19A mutant 
lacking all 24 effectors can still show biotrophic growth in plant 
tissue and thus behaves like S. reilianum (except for the systemic 
spread). Interestingly, two of the Tin4 orthologs of S. reilianum, 
sr 10075 and sr 10077, were recently shown to suppress apical 
dominance after maize infection (H. Ghareeb, F. Drechsler, C. 
Lofke, T. Teichmann and J. Schirawski, personal communication). 
The effect on apical dominance is a late phenotype observed about 
six weeks after infection of maize plants with S. reilianum, i.e. a time 
point not covered by our assays. To ascertain whether Tin4, 
Srl0075 and Srl0077 have conserved functions it would be 
interesting to test whether tin4 of U. maydis can complement the 
apical dominance phenotype of the respective S. reilianum mutants. 

Tinl-1 to Tin 1-5 is a group of weakly related U. maydis effectors, 
which could not be functionally separated because their individual 
effects on virulence were too small to be reliably detected. The 
transcriptome changes of plants infected with a mutant lacking all 
five related genes revealed specific, strong inductions of endochi- 
tinases, SA-binding proteins and the apoplastic peroxidase 
POX12. POX12 was recendy shown to be inhibited by the U. 
maydis effector Pepl leading to a suppression of the PAMP- 
triggered oxidative burst [16]. In addition, an NBS-LRR class 
disease resistance gene (Zm.3568.1) that could be involved in 
PAMP perception, showed transcriptional induction specifically 
after SG200Al9A-lb infections. Together, these changes indicate 
an enhanced defense response against the 19A-lb deletion mutant, 
which suggests that the Tinl-1 to Tin 1-5 effectors contribute to 
the suppression of basal host immunity. Interestingly, the immune 
response triggering avirulence factor, UhAvrl (UHOR_10()22, 
Figure 2), of U. hordei is most closely related to the U. maydis effector 
Tinl-2 and Tinl-3 [31]. With respect to virulence no specific 
contribution of UhAvrl could be detected [31], which may be 
consistent with the very small contribution to virulence that is seen 
when all five tinl genes are deleted. 

Tin3 is the effector in cluster 19A, which contributes most 
strongly to virulence. The strong transcriptional induction of two 
sucrose synthases after infection with the tin3 mutant strain (as well 
as after infection with the cluster 19A mutant) is likely to reflect 
enhanced photosynthetic activity in contrast to infections with 
SG200 where the transition from a juvenile sink tissue to a mature, 
photosynthetically active source tissue is blocked in infected leaves 
[9] . If an interplay between sucrose and auxin signaling, which was 
established in A. thaliana [32] also exists in maize; this could explain 
the observed upregulation of several auxin response factors after 
infection with these mutant strains. The specific upregulation of a 
WRKY transcription factor after infection with the cluster 19A 
mutant as well as the lin3 mutant could indicate elevated defense 
responses [33], which are downregulated by Tin3 after infections 
with wild type strains. Alternatively, this regulatory gene might 
negatively control cell cycle and/ or cell expansion, a feature of U. 
maydis induced tumors [9] . Another gene exclusively upregulated 
after infection with strains deleted for tin3 or cluster 19A is 
cytokinin oxidase 3, an enzyme involved in cytokinin degrada- 
tion. Cytokinin oxidases have been shown to restrict cell division 
and to regulate the sink capacity of kernels [34]. Thus, the 
downregulation of these activities by Tin3, presumably after 
uptake of Tin3 by plant cells, might be necessary for tumor 
development. 

The finding of discrete plant responses after infection with 
individual effector mutants provides important leads for the 
functional analyses that can now be followed. For example, the 
predicted changes in hormone levels attributed to Tin3 could be 
determined from metabolic profiles or directly connected with 
Tin3 by transiendy expressing Tin3 in plants with appropriate 



reporter gene constructs. The expression of Tin3 in transgenic 
plants might even allow assessing, whether the predicted effects on 
photosynthesis are direct or indirect. 

In more general terms our analyses reveal the power of studying 
pathogen effector mutants, by combining virulence assays with an 
assessment of plant responses to these mutants. Such comparisons 
do not only reveal common plant responses that reflect central 
processes targeted by the infection but in addition provide specific 
leads to the function of individual effectors. Furthermore, this 
approach does not rely on a significant virulence phenotype of the 
effector mutants studied and may thus be highly useful for the 
analysis of the vast majority of eukaryotic pathogen effectors that 
fall into this class [35,36]. 

Materials and Methods 

Fungal strains, growth conditions and plant infections 

U. maydis strains were grown in YEPSL (0.4% yeast extract, 
0.4% peptone, 2% sucrose) with shaking at 28°C at 200 rounds 
min 1 (rpm), to an optical density (OD 600 ) of 0.6-0.8. Cells were 
centrifuged at 3000 g for 5 min, resuspended in H z O to an OD 600 
of 1 and used for syringe infection of 7-day-old maize seedlings 
(variety Early Golden Bantam, Olds Seeds, Madison). At least 3 
biological replicates were tested for virulence and disease was 
scored 12 dpi following described protocols [12]. To demonstrate 
the statistical differences of disease symptoms in the mutants 
compared to SG200 each of the symptoms was tested by Student's 
Hest (** />-values<0.01.) and corrected by Bonferroni correction 
for multiple testing. The haploid solopathogenic strain SG200 [12] 
was used for virulence assays and all mutations were introduced 
into this background. 

Strain construction 

Standard molecular cloning strategies and techniques were 
applied [37]. All U. maydis strains (Table S5) are derived from the 
solopathogenic strain SG200 and were generated by a PCR-based 
gene replacement approach using primers listed in Table S6 or, for 
complementation experiments, by insertion of pi 23 derivatives 
into the ip locus as described [38]. Deletion endpoints are depicted 
in the respective Figures. Constructs used for complementation 
contained the respective gene plus the promoter region extending 
up to the next gene plus the Tnos terminator. All generated 
constructs were sequenced prior to U. maydis transformation (Table 
S7). Isolated U. maydis transformants were tested for integration 
events in the desired loci by southern analysis. For the 
complementation constructs, single copy integrations into the ip 
locus were selected by southern analysis. For 3 strains only 
derivatives containing two inserts could be obtained, this is marked 
in the strain list (Table S5). 

Quantitative real time PCR 

For quantification of relative fungal biomass in infected maize 
leaves 7-day-old maize seedlings were infected with SG200, 
SG200A19A, SG200A19A-1 and SG200A19A-2 and a section of 
the third leaf between 1 and 3 cm below the injection site was 
harvested after 2 dpi and 6 dpi. For genomic DNA extraction leaf 
material was frozen in liquid nitrogen, ground to powder, and 
extracted using a phenol-based protocol modified from Hoffman 
and Winston [39]. The qRT-PCR analysis was performed using 
an iCycler (Bio-Rad) in combination with the Platinum SYBR 
Green Supermix (Invitrogen). U. maydis biomass was quantified 
with primers PPI-fw (5 ' - AC ATC GTC AAGGC T ATC G- 3 ' ) and 
PPI-re (5 ' - AAAG AAC AC C GG AC TTGG- 3 ' ) amplifying the fun- 
gal ppi gene. Maize glyceraldehyde dehydrogenase was amplified 
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Figure 6. Graphic representation of Gene Ontology terms. Hierarchical presentation of Gene Ontology [43] terms showing molecular 
functions significantly enriched amongst upregulated maize transcripts in tin mutant infected tissue. Numbers give p-values for enrichment of the 
respective GO. p1 : p-value of enrichment in upregulated genes after infection by SG200Atin3 (pink boxes); p2: p-value of enrichment in upregulated 
maize genes after infection by SG200Atin4 (green box); p3: p-value of enrichment in upregulated maize genes after infection by SG200Atin5 (blue 
boxes). Box in yellow indicates significant enrichment in samples p1 and p2; box in dark pink indicates significant enrichment in p1 and p3; boxes in 
turquoise indicate significant enrichment in p2 and p3 and darker color symbolizes more significant enrichment. Red arrows indicate hierarchical 
connections between (significantly enriched) Gene Ontology terms. Dotted arrows represent GO terms not enriched and therefore spared out from 
the figure. 

doi:10.1371/journal.ppat.1003866.g006 



with primers GAPDH-F (5'- CTTCGGCATTGTTGAGG- 
GTTTG-3') and GAPDH-R (5'- TCCTTGGCTGAGG- 
GTCCGTC-3') [15] and served as reference gene for normali- 
zation. Relative amounts of fungal DNA (ppi) were then calculated 
relative to the amount of GAPDH DNA using the cycle threshold 
(Ct)2" AACt method [40]. Three biological replicates were com- 
bined and ^-values were determined by using Student's i-test (** p- 
values<0.01.). 

To validate the expression data of the microarray experiment 
13 maize genes differentially regulated after infection with 
different mutant strains were subsequently analyzed by 



qRT-PCR. Infected plant material was generated as described 
for the microarray experiment and used for RNA extraction with 
Trizol (Invitrogen, Karlsruhe, Germany). After extraction, the 
first-strand cDNA synthesis kit (Invitrogen) was used to reverse 
transcribe 3 |Xg of total RNA with oligo(dT) Primers. The qRT- 
PCR analysis was performed using an iCycler (Bio-Rad) in 
combination with the SYBR Green Supermix (Invitrogen). 
Primers used for quantification of maize gene transcription levels 
are listed in Table S6. Gene expression levels were calculated 
relative to GAPDH expression levels using the cycle threshold 
(Ct)2" AACt method [40]. 
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DNA microarray analysis 

For the microarray experiments, maize plants (Early Golden 
Bantam) grown under defined conditions in a growth chamber 
were infected with SG200A19A-lb, SG200Atin2, SG200Atin3, 
SG200Atin4 and SG200Atin5 as described previously [15]. 
Samples of infected tissue were collected 4 dpi by excising a 
section of the third leaf between 1 and 3 cm below the injection 
site. For RNA extraction, material from >20 plants per 
experiment was combined, ground to powder on constant liquid 
nitrogen and RNA was extracted with Trizol (Invitrogen, 
Karlsruhe, Germany). RNA was purified applying the RNeasy 
kit (Qiagen, Hilden, Germany). Affymetrix maize genome 
microarrays were performed in three biological replicates, using 
standard Affymetrix protocols (Midi_Euk2V3 protocol on Gene- 
Chip Fluidics Station 450; scanning on Affymetrix GSC3000G). 
Expression data were submitted to GeneExpressionOmnibus 
(http://www.ncbi.nlm.nih.gov/geo/) (Accession Number: 
GSE48406). 

Previously published Affymetrix data for SG200 infections [9] 
(GEO accession Number: GSE 10023) and the microarrays 
performed in this study were analyzed together using the Partek 
microarray software suite version 6.12. Expression values were 
normalized using the RMA method. Criteria for significance were 
a corrected /(-value (per sample) with a FDR of 0.05 and a fold- 
change of >2. Differentially expressed genes were calculated by a 
1-way ANOVA model using method of moments [41]. 

Microscopy 

For microscopic analysis of different life cycle stages of U. maydis 
strains SG200 and SG200A19A, a section of the third leaf between 
1 and 3 cm below the injection site was excised after 1 dpi, 13 dpi 
and 30 dpi. We used a Zeiss Axiophot with differential interfer- 
ence contrast (DIC) optics for microscopic observations. The 
pictures were taken using a CCD camera (C4742-95, Hama- 
matsu). To visualize penetration events, appressoria were stained 
with calcofluor white (100 |J,g/ml; Fluorescent Brightener 28, 
Sigma-Aldrich, Deisenhofen) for 1 min. Intracellular growing 
fungal hyphae were stained with chlorazol black E using an 
established protocol [42]. 

Supporting Information 

Figure SI Amino acid sequence alignment of gene families on 
cluster 19A. Amino acid sequences of gene families located on 
cluster 19A were aligned by Clustal Omega program. Similar 
amino acid residue among all was highlighted as black box. 
Similar amino acid residue among several amino acid sequences 
was highlighted by grey. (A) Orange group, (B) Yellow group, (C) 
Blue group, (D) Green group, (E) Brown group, described in 
Figure 1A. 
(PPTX) 

Figure S2 Quantitative real-time PCR of the genes upregulated 
in the leaves infected with tin mutants. The genes upregulated in 
the leaves infected with tin mutants were picked up from 
microarray data. RNA samples were extracted from the leaves 
inoculated with H z O or infected with SG200 and tin mutants, 
which were prepared independently, at 4 dpi, and used for qRT- 
PCR. Error bars were calculated from three biological replicates. 
(PPTX) 

Figure S3 Gene ontology enrichment analysis of maize genes 
induced by the cluster 19A deletion but not by the tin mutants at 
4 dpi. The GOEAST software toolkit [43] was used to identify 
GO terms for cellular processes (yellow boxes) that are specifically 



enriched in maize leaves infected with U. maydis strain 
SG200A19A samples. Darker color shades indicate higher 
significance of enrichment, /(-values are indicated in brackets. 
(PPTX) 

Figure S4 Gene ontology enrichment analysis of maize genes 
downregulated by the cluster 19A deletion but not by the tin 
mutants at 4 dpi. The GOEAST software toolkit [43] was used to 
identify GO terms for cellular processes (yellow boxes) that are 
enriched in maize leaves infected with U. maydis strain 
SG200A19A. Darker color shades indicate higher significance of 
enrichment, /(-values are indicated in brackets. 
(PPTX) 

Figure S5 Gene ontology enrichment analysis of maize genes 
downregulated by the tin mutants and the cluster 19A deletion at 
4 dpi. The GOEAST software toolkit [43] was used to identify 
GO terms for cellular processes that are enriched in maize leaves 
infected with U. maydis strain SG200A19A and the tin mutants 
compared to mock-treated samples. Yellow boxes indicate 
processes enriched in both SG200A19A and the tin mutants 
infected samples. Darker color shades indicate higher significance 
of enrichment, /(-values are indicated in brackets. 
(PPTX) 

Figure S6 Gene ontology enrichment analysis of maize genes 
induced after infection with U. maydis strain SG200Atin3 at 4 dpi. 
The GOEAST software toolkit [43] was used to identify GO terms 
for cellular processes (yellow boxes) that are specifically enriched in 
maize leaves infected with U. maydis strain SG200Atin3 Darker 
color shades indicate higher significance of enrichment, /(-values 
are indicated in brackets. 
(PPTX) 

Figure S7 Gene ontology enrichment analysis of maize genes 
induced after infection with U. maydis strain SG200Atin4 at 4 dpi. 
The GOEAST software toolkit [43] was used to identify GO terms 
for cellular processes (yellow boxes) that are specifically enriched in 
maize leaves infected with U. maydis strain SG200Atin4. Darker 
color shades indicate higher significance of enrichment, /(-values 
are indicated in brackets. 
(PPTX) 

Figure S8 Gene ontology enrichment analysis of maize genes 
induced after infection with U. maydis strain SG200Atin5 at 4 dpi. 
The GOEAST software toolkit [43] was used to identify GO terms 
for cellular processes (yellow boxes) that are specifically enriched in 
maize leaves infected with U. maydis strain SG200Atin5. Darker 
color shades indicate higher significance of enrichment, /(-values 
are indicated in brackets. 
(PPTX) 

Table SI Differentially regulated maize genes after infection 
with SG200, SG200A19A and SG200-derived tin mutant strains. 
(XLSX) 

Table S2 Differentially expressed maize genes shared and not 
shared in response to the cluster 1 9A deletion mutant and different 
tin gene deletion mutants. 
(XLSX) 

Table S3 Differentially regulated maize genes specific to 

infection with individual tin mutant strains. 

(XLSX) 

Table S4 List of Gene Ontology terms and enriched transcripts. 
(XLSX) 

Table S5 U. maydis strains used in this study. 
(DOCX) 
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Table S6 PGR primers used in this study. 
(DOCX) 

Table S7 Plasmids used in this study. 
(DOCX) 
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